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■ Abstract. I review here the history of bottom quark cross section measurements and theoretical 

^ \ predictions. Starting form the early days of UAl, and going through the sequence of the large 

Q ■ excesses reported during run and I at the Tevatron by CDF and DO, I summarize how both data and 

' theory have evolved in time, thanks to improved experimental techniques, more data, and improved 

control over the main ingredients of the theoretical calculations. I conclude with the discussion of 
the preliminary data from run II, which appear to finally give a satisfactory picture of the data vs 
theory comparison. 
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g : INTRODUCTION 
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The study of events with bottom quarks has led in the past 10 years to some of the most 
important Tevatron results: the discovery and study of the top quark, the appreciation 
of the colour-octet-mediated quarkonium production mechanisms, as well as general 
results in b-hadron physics (spectroscopy, lifetimes, mixing, sin^ 2/3) These results have 
, been obtained while both CDF and DO were reporting factor-of-3 discrepancies between 

^1 observed and predicted Z?-hadron cross-sections. To claim that we need to understand 

^ ■ b production in order to make new discoveries is therefore a bit exagerated: important 

discoveries should be able to stand on their feet without appealing to the prediction 
k>( ; of a QCD calculation. Nevertheless, lack of confidence in the ability to describe the 

j_] ■ properties of events containing b quarks, in addition to raising doubts over the general 



applicability of perturbative QCD in hadronic collisions, does limit our potential for the 
observation of new dynamical regimes (e.g. small-x physics |[ll]-||3]) or for the discovery 
of new phenomena (e.g. Supersymmetry [5]). In some cases, the existing measurements 
challenge the theory in ways which go beyond simple overall normalization issues, 
pointing at effects which are apparently well beyond reasonable theoretical systematics: 
this is the case of recent CDF studies, which detected anomalies in both rates and 
properties of events with secondary vertices and soft leptons |6]. It cannot be contested, 
therefore, that the study of b production properties should be one of the main priorities 
for Run II at the Tevatron, with implications which could go beyond the simple study of 
QCD. 

Starting from the situation as it developed during the early Tevatron runs, I will review 
here the progress in the theoretical predictions. More details on the historical evolution 
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of the cross section measurements can be found in fT], as well as in Jsll^, which also 
review the status of fixed- targ et heavy quark studies. For a recent review including 77 
and ep data as well, see 11 1011 . I will then present the implications of the preliminary 
results from Run II. Their complete theoretical analysis is contained in 



REVIEW OF RUN AND RUN I RESULTS 

The prehistory of b cross-section measurements in hadronic collisions starts with UAl 
at the SppSU/S = 630 GeV) collider [12]. The data were compared with theoretical 
predictions iflsL [l4ll . showing good agreement, within the rather large (±40%) theo- 
retical uncertainty. "Theory", in those days, already meant a full NLO QCD calcula- 
tion OEl], including all mass effects, state-of-art NLO PDF fits lUsll . and b-*B non- 
perturbative fragmentation functio ns p arameterized according to flu], with a parameter 
£ = 0.006 extrapolated from fits (ivlj to charm fragmentation data in e+e , using the 
relation = £c x (mc/mt,)^. At the beginning only predictions for total cross-sections 
and inclusive pj spectra were available. Later on, more exclusive calculations were per- 
formed, allowing for the application of general cuts to the final states, as well as for the 
study of correlations between the b and b [18]^. 

After such a good start in UAl, the first published data from CDF ifioll appeared as a 
big surprise. CDF collected a sample of 14 ±4 fully reconstructed XjrK^ decays, 

leading to: 

oipp m Pi > 11.5GeV. M < 1) = : 6,1 ± 1.9- ±2.V,Mb 

In spite of the large uncertainties, which led to a mere 1.5 a discrepancy, attention 
focused on the large data/theory =5. 5 excess. Theoretical work to explain the apparent 
contradiction between the success of the NLO theory at 630 GeV and the disaster at 
1.8 TeV concentrated at the beginning on possible effects induced by the different x 
range probed at the two energies: PDF uncertainties ll20ll and large small-jc effects 

where x ~ m\,l \fS. In the first case marginal fits to both data sets could be obtained 
at the cost of strongly modifying the gluon density, in a way which however would 
not survive the later accurate determinations of g{x) from HERA. In the second case, 
conflicting conclusions were reached: on one side the first paper of 10] obtained increases 
by factors of 3-5 due to small-x effects; on the other, the analysis of [T] proved that the 
resummation of small-jc logarithms could only augment the total rate by 30% (or less, in 
the case of g{x) more singular than 1 jx)^ . 

The ball was therefore back on the experimentalists' court. CDF expanded the set 
of measurements, including final states with inclusive i// and i//' [22] and inclusive 
leptons |23], summarised in fig. [U The measurement of the b cross section from the 



For lack of time, I will however focus my attention in this presentation on inclusive pj spectra. 
^ The option of very large small-jc effects being manifest only at 1 .8TeV will be definitely ruled out several 
years later, when CDF measured I21I1 the h cross section at \fS = 630GeV and showed that the scaUng 
from 630 to l.STeV was consistent with the predictions of pure NLO QCD. 
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FIGURE 1. CDF data from inclusive Xj/, y/' |22] and lepton |23] final states, compared to NLO QCD. 



inclusive charmonium decays turned out later to be incorrect. In run 0, in fact, CDF 
could not measure secondary vertices, so that charmonium states from direct production 
and from B decays could not be separated. The extraction of a b rate from these final 
states was based on theoretical prejudice about the prompt production rates, prejudice 
which in run I, when the secondary vertices started being measured by CDF, turned out 
to be terribly wrong [24]"^. The data on inclusive leptons, while high compared to the 
central value of the theoretical prediction, were nevertheless consistent with its upper 
value, and in any case within la. 

Increased statistics in run I allowed CDF to improve its measurement of fully recon- 
structed exclusive decay modes, leading to the measurements in fig. |2l For this mea- 
surement CDF used 19pb^^ of data, leading to approximately 55 — > x^rK* and 125 
B^ xirK^ decays. The cross section was still high compared to the central value of the 
theoretical prediction (data/theory=1.9±0.3), but this was already a marked improve- 
ment over the first measurement from run 0, when this ratio was equal to 6.1 ! More ex- 
plicitly, the 1995 measurement gave o{pt{B+) > 6GeV, \y\<\)= 2.39±0.54/ib, com- 
pared to the 1992 measurement of {o{pt{B) > 9GeV, \y\ < 1)) = 2.8 ± 1.4/ib (where 
{(y{B)) = [a (5+) +o{B'^)]/2). Taking into account that the rate is expected to increase 
by 2.7 when going from a 9 GeV to a 6 GeV threshold, the 1992 measurement appears 
to be a factor of 3.2 higher than the 1995 result, consistent with the 6.1/1.9 ratio. This 
drop in the experimental cross-section was not inconsistent with the large statistical and 
systematic uncertainties of the 1992 measurement, but somehow the common belief that 
theory was way off had already stuck. It is also worth noting that the same data, when 
compared to theoretical predictions obtained a couple of years later using the same QCD 
calculations, but up-to-date sets of input PDFs (MRST [27] with aj(mz)=0.1175, and 
CTEQ5M ll28l] with as{mz)=0.1lS), gave very good agreement. This is shown in the 



Incidentally, this fact puts into question the UAl results, which heavily relied on the \j/ final states and 
on explicit assumptions about the prompt charmonium rates! 
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FIGURE 2. Evolution of data/theory comparisons with improved PDF fits. The data on both plots are 
exactly the same; the theory curves on the left were generated with the MRSDO set, on the right with the 
post-HERA set CTEQ5 and MRST. 
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FIGURE 3. Left: the NLO Zj-quark rate as a function of pr.min, for post-HERA PDF sets CTEQ4M 
(|29|, a,(mz)=0.116) and CTEQ6M (|301, a,(mz)=0.1 18), normaHzed to the pre-HERA set MRSDO. 
Right: total cross section for \y\ < 1 for various PDF sets, distributed on the abscissa in order of increasing 
release date. The crosses correspond to the rates calculated by forcing Aqcd to take a value consistent with 
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the LEP a,(mz) fits (A^-i:^ : 



: 226MeV => a,(mz) = 0.118 



right panel of fig. |2l taken from an update of ||^. The crucial change between the two 
predictions was the change in the value of the QCD coupling strength tts extracted from 
global PDF fits. The fits used in the CDF 1995 pubUcation, MRSDO |26], did not include 
HERA data and had as{mz)=0.lll, signficantly lower than what we were getting from 
LEP, namely as(mz) ~ 0.120. This 10% difference, when evolved to the low scales of 
relevance to b production, becomes much more significant, especially because b rates 
grow like aj. This is shown more explicitly in fig.|3 The left panel shows the ratio of 
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FIGURE 4. Final CDF analysis of run I exclusive-decay data fsH, compared to the CDF evaluation of 
the NLO QCD prediction with MRST PDFs and Peterson fragmentation. 



the rates obtained by using post-HERA PDF sets, normalized to the MRSDO set used 
in the CDF 1995 comparison. The right panel shows the integrated total cross section 
(for \y\ < 1) for several PDF sets, ordered versus the date of release. One can notice a 
constant increase, with the most recent sets being almost a factor of 2 higher than the 
older ones. Notice that this increase is due by and large to the increased value of tts 
returned by the PDF fits. Forcing Aqcd to take the value consistent with LEP's as{mz), 
one would have obtained for each PDF set the values corresponding to the crosses in the 
plots. There the increase relative to the pre-HERA fit MRSDO is significantly smaller. 

While the improvements in the PDF fits were reducing the difference between data 
and theory, as shown fig. |2l a new CDF measurement from the full sample of run I 
exclusive B decays in the range 6 GeV < pr <20 GeV appeared in 2001 llsill . and is 
shown here in Fig.|4| The total rate turned out to be 50% larger than in the previous 1995 
publication 125 1: o{pt{B^) > 6GeV, \y\ < 1) = 3.6±0.6/ib, compared to the previous 
2.4 ±0.5jUb, a change in excess of 2a. The ratio between data and the central value of 
the theory prediction was quoted as 2.9 ±0.5: a serious disagreement was back! 

On the other side of the Tevatron ring, the DO experiment started presenting the 
first b cross section measurements in 1994. The first preliminary results [32] were in 
perfect agreement with QCD, as shown in the left panel of Fig.|5l They were eventually 
published, after significant changes, in 113 3h . The results from a larger dataset of 6.6pb^^ 
appeared in [34], where i// dimuons were added. They are shown in the central panel 
of the figure, and they show a clear increase over the preliminary analysis, but are still 
consistent with the QCD expectations. The same data set underwent further analysis, 
and eventually appeared few years later in [35]. They are shown in the right panel of 
the figure. Now the data are significantly higher than QCD, and certainly higher than 
in 1996, especially in view of the fact that in the meantime the theory predictions 
had increased by almost a factor of 2 as a result of the use of new PDF sets (this is 
clearly visible by the shift of the theory curves between the central and right panels). 
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FIGURE 5. Evolution of the DO measurements. Left: preliminary results from 90nb-i HI. Center: 
6.6pb^^ |34|. Right: final analysis of the same data set, with the addition of inclusive dimuons f35ll. 
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FIGURE 6. Forward muon production at DO (3; 




As in the case of the CDF exclusive analysis, this evolution underscores the difficulty 
in performing these measurements, and indicates that it was not just the theory that was 
having difficulties in coming to grips with the problem! 

An additional element was added to the puzzle when DO reported ll36ll the measure- 
ment of b production at large rapidity, using inclusive forward muons (2.4 < |j;x | < 3.2). 
The results, shown in fig. |6l indicated an excess over NLO QCD by a factor larger 
than what observed in the central region. This anomaly could not be explained away by 
assuming some extra systematics related to PDFs. From the point of view of perturba- 
tion theory, furthermore, there was no reason to expect a significant deterioration of the 
predictive power when going to large rapidity. So when this result first appeared in its 
preliminary form I was led [38] to review our assumptions about the non-perturbative 
part of the calculation, in particular the impact of the fragmentation function. A crucial 
observation is that in hadronic collisions the fragmentation function is probed in differ- 
ent ranges of z as we change rapidity. This is easily seen as follows. Let us assume that 
the b pt spectrum takes the simplified form: 



doib) 
dpT 



1 



(2) 



where the slope A'^ will typically depend on rapidity, becoming larger at higher j^. The 
meson spectrum is then obtained via convolution with the fragmentation function f{z). 




20 JO 40 50 60 70 BO 90 100 



E,'''' (GeV) 

FIGURE 7. fe-jet production at DO Ii40ll. 



leading to the simple result: 

do{B) _ [dzdo{b) [^^(^\Nri^ r do{b) 

-\Pt=Pt/z)= —[—) f{z)=fN ,„ , (3) 



dPr J z dpT J z Pt dPr 

where /n is the N-th moment of f{z) ■ This means that a steeper partonic spectrum selects 
higher moments. Since the index A'^ is larger for forward production, a relative difference 
in B production rates in the forward/central regions could be explained by making the 
fragmentation function harder, enhancing the larger moments (which measure the large- 
z behaviour of f{z)). A related observation is that f(z) fits to e^e^ data are mostly driven 
by the value of the first moment /i , which measures the average of the fragmentation 
variable z. It is therefore possible that different choices of /(z), giving equivalent overall 
fits to e+e^, could make very different predictions for the higher moments of relevance 
to hadronic production (in this case A'^ is in the range 4-6). 

One way to understand whether indeed the inaccurate description of the fragmentation 
process could affect the theoretical predictions was therefore to think of measurements 
not affected by this systematics. The most obvious observable of this kind is the Et 
spectrum of jets containing a b quark Since the tagging of ab inside the jet is only 
marginally affected by the details of the b ^ B fragmentation, measuring the rate of 
b jets is a direct measurement of the b production rate with negligible fragmentation 
systematics. In addition, this measurement is also insensitive to higher-order large- 
Pt logarithms which are present in the pj spectrum, therefore improving in principle 
the perturbative accuracy. DO carried out the measurement, publishing [40] the results 
shown in Fig. IT] The agreement with NLO QCD [39] is better than in the case of the 
p^ spectrum, as was hoped. We took this as strong evidence that a reappraisal of the 
fragmentation function systematics may have led to a better description of the pj and 

distributions. 

The necessary ingredients to carry out this programme are perturbative calculations of 
matching accuracy for b spectra in both e+e^ and pp collisions, in addition of course to 
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FIGURE 8. Left panel: FONLL prediction by Cacciari and Nason 14511 for the run I B meson spectrum, 
compared to the CDF data |31]. Right panel: the prediction of this calculation for the forward muon 
rapidity spectrum at DO. 



accurate e~^e~ data to be used in the fits. These tools had just become available towards 
the end of the 90's The resummation of the logarithms of pr/mij, with next-to-leading 
logarithmic accuracy (NLL), and the matching with the fixed-order (FO), exact NLO 
calculation for massive quarks, had been performed in 114 ih (Fixed-Order with Next- 
to-Leading-Log resummation: FONLL) and a calculation with this level of accuracy 
for e+e^ collisions was presented in [42]. Here it had been used for the extraction of 
the non-perturbative fragmentation function f{z) from LEP and SLC data ll43ll . with 
the main result that the Peterson functional form is strongly disfavoured over other 
alternatives [44]. The equivalence of the perturbative inputs allows one to consistently 
apply this fit to the FONLL ^-quark spectra in hadronic collisions, leading to FONLL 
predictions for the b hadron (Hj,) spectrum. A comparison of these predictions with the 
final CDF data at 1 .8 TeV for 5^-meson production in the range 6 GeV < pr <20 GeV 
has been presented in 11451] . The results are shown in Fig. [HI the left panel compares the 
CDF data from fsi'] with the theory curve evaluated using CTEQ5M PDF, FONLL, 
and fragmentation functions fitted to LEP and SLC data. The right panel shows a 
comparison ll4^ with the DO forward muon data. In both cases the agreement with data 
is much improved. In the case of the CDF central cross section, the ratio between data 
and theory improves from 2.9 ±0.5 to 1.7 ±0.7. As discussed in detail in [45], the 
improvement is due to the sum of three independent 20% effects (1.23 _ 2.9/1.7), all 
going in the same direction: the resummation of pr logarithms, the change in functional 
form of the fragmentation function, and the use of the LEP/SLC b fragmentation data. 
The heritage of run I was therefore a set of measurements, more or less consistent with 
each other, normalized with a factor of about 1 .5 to 2 higher than the central theoretical 
prediction, but still compatible with the upper end of the theoretical systematics band. 
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FIGURE 9. CDF 7/^/ spectrum from B decays. The theory band represents the FONLL systematic 
uncertainties, as described in the text. Two MC@NLO predictions are also shown (histograms). 



THE RUN II CDF RESULTS 



The final phase of this history deals with the new run II data from CDF 114 811 . A great 
improvement took place in the ability to trigger on very low pj events, allowing for 
a measurement down to ~ 0, although still in the limited rapidity range \yi,\ ^ 0.6. 
This is also accompanied by very large statistics, allowing a fine binning in pj. The 
measurement down to very small pj is important because the total rate has a much 
reduced dependence on the fragmentation systematics, and because it is particularly 
sensitive to possible small-jc phenomena. 

On the theoretical side, in addition to the calculations described above, a new tool 
has meanwhile become available, namely the MC@NLO code ll4^ . which merges 
the full NLO matrix elements with the complete shower evolution and hadronization 
performed by the HERWIG Monte Carlo. As discussed in detail in [49], this comparison 
probes a few features where FONLL and MC@NLO differ by effects beyond NLO: the 
evaluation of subleading logarithms in higher-order emissions, in particular in the case 
of gluon emission from the b quark, and the hadronization of the heavy quark, which in 
MC@NLO is performed through HERWIG's cluster model, tuned on Hj^X decays. 

The comparison of the run II data with the theoretical calculations is given in Fig. |9l 
which shows the data with our prediction for the spectrum of 7/i//'s form Hi, decays, 
obtained by convoluting the FONLL result with the 7/i// momentum distribution in in- 
clusive B ^ J /x^r + X decays. The theoretical error band is obtained by varying renor- 
malization and factorization scales {IIr,f = ^r.fIM), with /Iq = pr^ + m^), the ^-quark 
mass, and parton densities. The central values of our predictions are obtained with 
^R,F = ^, mb = 4.75 GeV and CTEQ6M. The mass uncertainty corresponds to the range 
4.5 GeV < < 5 GeV. The scale uncertainty is obtained by varying ji^ f over the range 
0.5 < ^R^f < 2, with the constraint 0.5 < ^r/^f < 2. The PDF uncertainty is calculated 
by using all the three sets of PDFs with errors given by the CTEQ, MRST and Alekhin 



groups 

The data lie well within the uncertainty band, and are in very good agreement with the 
central FONLL prediction. I also show the two MC@NLO predictions corresponding to 
the two different choices of the b hadronization parameters (see iflll] for the details). 

I stress that both FONLL and MC@NLO are based on the NLO result of [fli 
(henceforth referred to as NDE), and only marginally enhance the cross section predicted 
there, via some higher-order effects. The most relevant change in FONLL with respect 
to old predictions lies at the non-perturbative level, i.e. in the treatment of the b 
Hh hadronization, which makes use [45] of the moment-space analysis of the most 
up-to-date data on b fragmentation in e+e annihilation. The evolution of the NLO 
theoretical predictions over time is shown in Fig.[lOl Here we plot the original central 
prediction of NDE for =1.8 TeV (symbols), obtained using NLO QCD partonic 
cross sections convoluted with the PDF set available at the time, namely DFLM260 [ fsl]- 
The same calculation, performed with the CTEQ6M PDF set (dotted curve), shows 
an increase of roughly 20% in rate in the region pr < 10 GeV. The effect of the 
inclusion of the resummation of NLL logarithms is displayed by the dashed curve, 
and is seen to be modest in the range of interest. Finally, we compare the original 
NDE prediction after convolution with the Peterson fragmentation function (e = 0.006, 
dot-dashed curve), with the FONLL curve convoluted with the fragmentation function 
extracted in ll45h (solid curve). Notice that the effect of the fragmentation obtained 
in 1I45I1 brings about a modest decrease of the cross section (the difference between 
the dashed and solid curves), while the traditional Peterson fragmentation with e = 
0.006 has a rather pronounced effect (the difference between the symbols and the dot- 
dashed curve). Thus, the dominant change in the theoretical prediction for heavy flavour 
production from the original NDE calculation up to now appears to be the consequence 
of more precise experimental inputs to the bottom fragmentation function [43], that 
have shown that non-perturbative fragmentation effects in bottom production are much 
smaller than previously thought. 

The main improvement in the comparison between data and theory w.r.t. the final 
run I results discussed in comes from the normalization of the run II CDF data, 
which tend to be lower than one would have extrapolated from the latest measurements 
at 1.8 TeV. To clarify this point, we collect in Fig. [TT] the experimental results from the 
CDF measurements of the cross section in Run lA lEsll . in Run IB ll3lll and in Run II. 
The rate for pt{B^) > 6 GeV, evolved from 2.4 ±0.5 (Run lA) to 3.6 ±0.6 
(Run IB), and decreased to 2.8 ±0.4 jih in Run II. The increase in the cm. energy 
should have instead led to an increase by 10-15%. The Run II result is therefore lower 
than the extrapolation from Run IB by approximately 30%. By itself, this result alone 
would reduce the factor of 1.7 quoted in [45] to 1.2 at =1.96 TeV. In addition, the 
results presented in 111] lead to an increase in rate relative to the calculation of [45] 
by approximately 10-15%, due to the change of PDF from CTEQ5M to CTEQ6M. 
We then conclude that the improved agreement between the Run II measurements and 
perturbative QCD is mostly a consequence of improved experimental inputs (which 
include up-to-date and PDF determinations). 
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FIGURE 10. Evolution of the NLO QCD predictions over time, for \/5 = 1800 GeV. See the text for 
the meaning of the various curves. 
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FIGURE 11. Evolution of the CDF data for exclusive B± production: Run IAf2?l, Run IB fs?) and 
Run IIL4SJ. 



CONCLUSIONS 

When I meet colleagues and discuss the latest b results, and when I hear presentations 
or read conference proceedings, there is often a more or less explicit message that now 
things are OK because theorists kept beating on their calculations until they got them 
right. I hope that this note will dispel this prejudice. The history of the experimental 
measurements indicates that many things have also "strongly evolved" on the data side, 
often with changes well in excess of the standard ±la variation. The "history" plot 
in fig. [jni shows on the other hand that not much has changed on the theory side, 
aside from data-driven modifications associated to the value of as{mz), to the low-x 



behaviour of the gluon as determined by the HERA data, and to the improved data 
on — i> 5 fragmentation. The theoretical improvements due to the resummation of the 
\arge-pT logarithms play a major role in allowing a consistent use of the fragmentation 
functions extracted from e+e data, but have a very limited impact in the region of 
probed by the run II data. Their significance will only manifest itself directly at high 
{pj > 20 — 30GeV), where the resummation leads to a much reduced scale dependence, 
and to more accurate predictions, allowing more compelling quantitative tests of the 
theory. It is auspicable that the improved run II detectors and the higher statistics will 
make it possible to extend the range of the measurements to really large pj (in the 
range of 80-100 GeV). Tools are now available (MC@NLO) to compare data subject 
to complex experimental constraints directly with realistic NLO calculations, including 
a complete description of the hadronic final state. This will avoid the risky business of 
attempting to connect the observables to a pt spectrum of the b quark, a practice which, 
although unavoidable in the past, has certainly contributed to the inflation of theoretical 
and experimental systematic uncertainties. 

To this date, the recent CDF measurement of total Z^-hadron production rates in 
pp collisions at \/S = 1.96 TeV is in good agreement with NLO QCD, the residual 
discrepancies being well within the uncertainties due to the choice of scales and, to a 
lesser extent, of mass and PDF. A similar conclusion is reached for the pj spectrum. 
The improvement in the quality of the agreement between data and theory relative to 
previous studies is the result of several small effects, ranging from a better knowledge of 
fragmentation and structure functions and of a^, which constantly increased in the DIS 
fits over the years, to the fact that these data appear to lead to cross sections slightly lower 
than one would have extrapolated from the measurements at 1.8 TeV. The currently 
still large uncertainties in data and theory leave room for new physics. However there 
is no evidence now that their presence is required for the description of the data, and 
furthermore the recent results of |52] rule out the existence of a scalar bottom quark in 
the range preferred by the mechanism proposed in The data disfavour the presence 
of small-.x: effects of the size obtained with the approaches of refs. iJ]. They are instead 
compatible with the estimates of [\S ■ 

While these results have no direct impact on other anomalies reported by CDF in 
the internal structure and correlations of heavy-flavoured jets [6], we do expect that 
the improvements relative to pure parton-level calculations present in the MC@NLO 
should provide a firmer benchmark for future studies of the global final- state stucture of 
bb events. 
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